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Abstract: The role of the Coulomb potential on the generation of ellip- 
tically polarized high-order harmonics from atoms driven by elliptically 
polarized laser is investigated analytically. It is found that the Coulomb 
effect makes a contribution to the harmonic ellipticity for low harmonic 
orders and short quantum path. By using the strong-field eikonal-Volkov 
approximation, we analyze the influence of the Coulomb potential on 
the dynamics of the continuum state and find that the obtained harmonic 
ellipticity in our simulation originates from the break of symmetry of the 
continuum state. 
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1. Introduction 

When atoms and molecules are exposed to intense laser irradiation, high-order harmonics (HH) 
will be generated. The high-order harmonic generation (HHG) has been an attractive topic for 
the past two decades for its potential application of producing coherent attosecond pulses (l}|7) 
in the XUV regime and self-probing of molecules allowing a combination of attosecond and 



Angstrom resolutions ||8]jl2) . The physical mechanism of HHG can be understood by the three- 
step model: (i) an atom or molecule emits its electron to continuous state, (ii) the electron is 
accelerated in the laser field, (iii) the accelerated electron recombines with the parent ion and 



high energy photon is radiated 1 13 14 1. As described by this model, linearly polarized high- 
order harmonics will be generated from atomic gases by linearly polarized laser pulses. It is 
also found that, when the laser pulses are elliptically polarized, elliptically polarized HH will 
be obtained fl5). In the observed phenomena, the rotation of the polarization of the HH field 



with respect to that of the driving laser can be explained by the classical three-step model 1 16 1, 
but this classical model can not explain the harmonic ellipticity. Quantum models based on 
the strong-field approximation (SFA) can be employed to provide theoretical approaches to 
investigate the observed ellipticity fT5][T7l[T8| . However, these theoretical approaches did not 
suggest qualitative explanation for the origin of the ellipticity. To answer the question what 
results in the harmonic ellipticity not only enables people to gain a deeper insight into the 
underlying physical mechanism of atomic HHG, but also help researchers understand more 
clearly about the generation of elliptically and circularly polarized molecular HH 1 19-24). 



In a recent work, the polarization properties of the HH generated by atomic gases in el- 
liptically polarized lasers was investigated and the ellipticity of the HH was attributed to the 
quantum-mechanical uncertainty of the electron motion J25| . The authors have developed a 
theoretical model to describe the ellipticity and their analytical result agrees well with the exact 
numerical result for HH of high orders, but the analytical result is a little lower than the numer- 
ical one for low harmonic orders (see Fig. 2 in E5J). In that work, the Coulomb effect is not 
taken into account and this effect is considered to be not responsible for the harmonic ellipticity 
at least for the short quantum path and not very low harmonic orders. In ref. [22 26 27 1, the 
authors have pointed out that the molecular potential is responsible for the observed ellipticity 
in molecular HHG, which inspires one that the deviation observed in ref. (25} for low harmonic 
orders may result from the neglect of the Coulomb potential. In ref. ]18[, the Coulomb effect 
was also introduced by using the continuum Coulomb state in the theoretical model for HHG in 
elliptically polarized field. However, in the previous works, how the potential affects the HHG 
and leads to the harmonic ellipticity is still not clear enough. Therefore, more detailed studies 
concerning the role of Coulomb effect on the generation of elliptically polarized HH are in 
demand. 

In this paper, the effect of the Coulomb potential on the generation of elliptically polarized 
HH from atoms driven by intense elliptically polarized lasers is investigated analytically. By 
using the strong-field eikonal-Volkov approximation (SF-EVA), we in detail analyze the influ- 
ence of the Coulomb potential on the dynamics of the continuum states and exposit the role 
of the Coulomb potential on the harmonic ellipticity. We further investigate the contribution 
of the Coulomb effect to the harmonic ellipticity for different harmonic orders and quantum 
paths, which explains the fact that the ellipticity is mainly obtained in the low order harmonics 
from short quantum path in our simulation. Moreover, the conclusions are further confirmed by 
comparing the results obtained with the soft-core potential and short-range potential. 

2. Theoretical model 

We employ 1300 nm, 2.2x 10 14 W/cm 2 , elliptically polarized driving lasers with ellipticity of 
0. 1 . The major axis of the polarization ellipse is parallel to the x axis of the laboratory frame and 
the minor axis is parallel to the y axis. To obtain the ellipticity of the HH, the time-dependent 
transition dipoles of the system are calculated by j2"8||29") 



di(t;k) = (\i/ (r,t)\ri\y k (r,t)),i = x,y 



(1) 



where y/b is the ground state of the atom and y/<. refers to the continuum state, k is the recom- 
bination momentum, which is associated with the recombination time t . 

In our simulation, the ground state y/o is obtained by solving the time-dependent Schrodinger 
equation (TDSE) with imaginary-time propagation. The Hamiltonian is (in atomic units) 

// = -V + V(r), (2) 

with 

V r z + a 

where the soft-core parameter a is set to be 0.388 to fit the ground-state energy of Ar E g = 
—0.5794 a.u. Through this paper, the target atom is Ar. 

Regarding the continuum state y/^, if the continuum electron is treated as a free particle 
moving in the electric field without the effect of Coulomb potential, the continuum state is 
approximated by the plane wave y/fc(r,£) = e' k M r . In the expression, the space-independent 
terms for the continuum state aj on [tb(t)]exp[— iJl tc^(t')/2dt']/(2n)^' 2 ' are omitted, where tj, 
is the ionization time, ki(f') is the instantaneous momentum of the continuum electron in the 
laser field and fl, on (^) refers to the ionization amplitude. These terms can be omitted because 
they will be canceled in calculating the amplitude ratio R and the phase different 8 between the 
two orthogonal components. The physical meaning of this treatment is that, according to the 
three step model, the HHG in the x,y components share identical ionization and acceleration 
steps and the different properties between the two components of the HH radiation originate 
from the recombination step 1 10 19 29) . 



For the plane waves, according to the Euler's formula e' k r — cos(k ■ r) + i sin(k ■ r), the real 
and imaginary parts of y/j are symmetric and antisymmetric on the spatial coordinates respec- 
tively. As the ground state y/Q is also symmetric, when its space-independent term e~' £ «' is 
also omitted, the result of the integral dj(t\k) — J y/or/y/^c/r must be a pure imaginary num- 
ber and the phase difference between the two components d y and d x is either or % |23| . In 
the more realistic situation, when the space-independent terms are included, the result for the 
phase difference remains the same. So the calculated HH is linearly polarized and the harmonic 
ellipticity cannot be explained using the plane wave approximation due to the symmetry of the 
states j2"9"l . 

In this work, we employ the SF-EVA 1 30- 32], which enables one to describe the continuum 
dynamics taking into account the Coulomb effect. With this approach, the continuum state is 
expressed as 

y/ fe (r,0=e ikW ' r+, ' <T * (r '' 3 . (4) 

a]t(r,f) is the correction to the plane wave considering the effect of the Coulomb potential and 
is given by 



a k (r,t) = -l k L (t'){k L (t')-^kl(t')-2V[r L (t')]}dt' (5) 

where 



r L (t')=T-J t k L (t")dt" (6) 

describes the motion of the continuum electron in the laser field. 

Compared with the continuum Coulomb state used in fl8) , which only depends on the mo- 
mentum k at time t, SF-EVA describes the Coulomb correction to the continuum state by trac- 
ing the motion of the electron in the process from % to t. Therefore, with this approach, it is 
possible to in detail investigate influence of the Coulomb potential on the dynamics of the con- 
tinuum electrons with different trajectories and to study the role of the Coulomb potential on 



the harmonic ellipticity. In ref. 1 3 1 1, a detailed and systematic introduction to the SF-EVA is 
presented. 

In this work, in the case of HHG driven by elliptically polarized lasers, we apply the same de- 
scription for the motion of the continuum electron as in (25]. In this description, the continuum 
electron starts from the origin with zero velocity and the return of the electron is judged by x=0. 
Therefore, the relation among t, % and k can be obtained by finding the approximate solutions to 
the saddle-point equations [ 17 33 - 35] in the x component. Then, the instantaneous momentum 
k/,(f) and the classical trajectory of the continuum electron can be calculated by substituting 
the solutions into the equations of the electron motion for both x and y components. 

The amplitude ratio and the phase difference between the two orthogonal harmonic compo- 
nents are calculated by R = \d y \/\d x \ and 8 — arg[d y ] — arg[d x ] respectively, and the ellipticity 
e is finally obtained by p6| 



/ 1 + R 2 - Vl + 2R 2 cos28 + R 4 

E = \ / (7) 

V l+R 2 + Vl+2R 2 cos28+R 4 

When the amplitude ratio R = or the phase difference 8 — or %, the ellipticity £ = and the 
HH is linearly polarized, otherwise it will be elliptically or circularly polarized. In this work, in 
order to study the role of Coulomb effect on the harmonic ellipticity, we only take into account 
the influence of the Coulomb potential in the model and the other effects such as the quantum 
mechanical uncertainty of the electron motion are not included. 



3. Result and discussion 

In Fig. 1, three typical trajectories of the continuum electrons with final recombination kinetic 
energy of IU P , 2U P and 3U P are presented, where U p — Fq/4-co 2 with Fq and (0 being the 
amplitude and angular frequency of laser. Panels (a) and (b) correspond to the short and long 
quantum paths respectively. The continuum electron starts from the origin with zero velocity 
and recombines to the core at x = 0, which is the same as in ref. (25). As discussed above, 
these trajectories are obtained by solving the saddle-point equations in the x component and 
substituting the approximate solutions into the equations of the electron motion for both x and 
y components. 

According to the cutoff law of HHG (33) , the recombination kinetic energy of 3U p corre- 
sponds to HH radiation at the cutoff in the harmonic spectrum, and the recombination kinetic 
energy of 2U P , IU P correspond to HH of lower orders in the plateau region (the 89th and 53rd 
orders respectively). As shown in Fig. 1(a) for short quantum path, when the recombination 
kinetic energy grows higher, the trajectory grows longer and the motion of the continuum elec- 
tron is farther from the core. The maximum distance increases from no more than 20 a.u. to 
nearly 60 a.u. As shown in Fig. 1(b) for long quantum path, the trajectories are even longer for 
HH from the cutoff to lower orders. At the same time, with the trajectory extending longer, the 
transverse shift grows up to 40 a.u. Figs. 1(c) and (d) are the snap shots of the evolution of the 
probability density distribution of the wavepacket at the instant of recombination. The evolu- 
tion is obtained by solving the TDSE numerically. In the quantum view of the dynamic of the 
continuum electron, one can see that although the center of the continuum wavepacket shifts 
about 30 a.u. transversely, it is still wide enough to interfere with the bound state to generate the 
HH. For the numerical evolution, the laser parameters are the same as those given in Section 
2 and the initial state is also obtained by solving TDSE with imaginary-time propagation with 
the Hamiltonian given by Eq. (2). 

It has been shown in Fig. 1 that the behaviors of the continuum electrons are different among 
different harmonic orders and quantum paths. These different features of the motions of the 
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Fig. 1. (a) and (b) Typical trajectories of the continuum electrons with recombination ki- 
netic energy of IU P , 2U P and 3U P for short and long quantum paths respectively, (c) and 
(d) Snap shots of the evolution of the probability density distribution of the wavepacket at 
the instant of recombination. 



continuum electrons will result in the different harmonic ellipticities in the HHG for different 
harmonic orders and quantum paths. According to the previous discussion, if the continuum 
state is approximated by the plane wave, the obtained HH will be linearly polarized. However, 
when we take into account the Coulomb effect, during the acceleration step of HHG, the contin- 
uum state will be modified by the Coulomb potential before recombination. This modification 
is expressed by the term a>(r,f) in Eq. (5), which is mainly determined by the integral for the 
Coulomb potential along the trajectories, i.e. V[ri(t')]. After considering this modification in 
the phase, the continuum state will be no longer a plane wave. As a result, the symmetry of the 
wave function % is broken and the obtained d x , d y are not pure imaginary numbers any more. 
Then, the value of the phase difference 8 varies in the interval of [0,27c) rather than fixes to 



or 7E 1 23 1 . Therefore, the Coulomb effect on the continuum states finally leads to the elliptically 
polarized HHG. 

In Fig. 2, we calculate the modification term (^.(r,?) at the instant of recombination with 
kinetic energy of IU P and 3U P and for short and long quantum paths respectively. All the four 
panels are plotted in the same color scale. Comparing the second column to the first column 
in Fig. 2, the Coulomb effect is weaker for long quantum path. This is because for the long 
quantum path, the electron leaves directly far away from the core, drifts in the faraway region 
where the Coulomb potential V(r) is weak and enters the vicinity of the core shortly before 
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Fig. 2. Calculations of the Coulomb modification term Cfy(r,/) at the instant of recombi- 
nation for (a) short quantum path and low harmonic order, (b) long quantum path and low 
harmonic order, (c) short quantum path and high harmonic order and (d) long quantum path 
and high harmonic order. 

recombination. While for short quantum path, the electron only drifts in the vicinity of the 
core, where the effect of the Coulomb potential is much stronger. 

In addition, comparing the first row to the second row, for both the short and long quantum 
paths, the Coulomb effect for lower harmonic orders is stronger than that for higher orders. 
This is because, in the case of short quantum path, the trajectory of the continuum electron 
for the lower harmonic order is closer to the core and the influence of the Coulomb potential 
V [ri(t )] is much bigger. Thus it is obvious that the Coulomb effect is strongest for the low order 
harmonics from short quantum path as shown in panel (a). For the long quantum path, although 
the trajectory for \U P is longer and farther from the core than that for 3U P , the Coulomb effect 
is already very weak for both faraway electrons. The determining factor is that, the accumulated 
Coulomb phase is greater for the slow electron especially in the processes when the electron 
leaves and comes back to the vicinity of the core. Mathematically speaking, the integral between 
tj, and t in Eq. (5) is bigger for HH of low harmonic order. 

Finally the calculated rotation angle of the HH field and the ellipticity as functions of the 
recombination time t are presented in Fig 3. The recombination at about t = 0.7 optical cycle 
is associated with the HH at the cutoff region in the harmonic spectrum, and recombination 
at earlier or later time corresponds to HH of lower harmonic orders in the plateau region. The 



rotation angle is defined by the direction of the electron momentum at the instant of return |25 1. 
Nearly identical curve is observed as that in ref. p5| , which indicates that our description to 
the motion of the continuum electron is the same. Regarding the ellipticity, as shown by the 
solid blue curve, unnegligible harmonic ellipticity is observed. The ellipticity is up to more 
than 0.2 for HH of low harmonic orders contributed by short quantum path. The result shows 
that the Coulomb effect also makes a contribution to the harmonic ellipticity at least for the low 
harmonic orders and short quantum path. 
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Fig. 3. The harmonic ellipticity and the rotation angle as functions of the recombination 
time /. 

When we transform the ellipticity curve as a function of the recombination time t into curves 
as a function of the harmonic orders, it will be more obvious for the result to be analyzed. This 
is shown in Fig 4(a) and the phase differences between the two orthogonal components 8 are 
also plotted in Fig 4(b). As the relationship between t and k has been obtained by solving the 
saddle-point equations in the x component, the recombination momentum is then transformed 
into the harmonic order q according to the relation q(0 = k 2 /2 — E g . 

It is found in Fig. 4(a) that the harmonic ellipticity of short quantum path is much higher than 
that of long quantum path, and that the ellipticity for both quantum paths drops as the harmonic 
order increases. According to the previous discussions and analyses, the obtained ellipticity 
originates from the break of symmetry of the continuum state and the break of symmetry re- 
sults from the Coulomb modification. In Fig. 2, it is shown that the Coulomb modification for 
the short quantum path is much stronger than that for the long quantum path, and the Coulomb 
modification for low harmonic orders is stronger than that for high harmonic orders. Therefore, 
all the observations above are consistent with the previous discussions and can be well ex- 
plained by the analyses of Fig. 2. The vanishment of the ellipticity for the high harmonic orders 
implies that the effect of the Coulomb potential on high energy electrons is much weaker and 
the Coulomb continuum state is reduced to approximate plane waves, which also agrees with 



the prediction in |25) . 

Comparing Fig. 4(b) with Fig. 4(a), the relationship between the ellipticity and the phase dif- 
ference is also consistent with the previous discussion. For high harmonic orders, the ellipticity 
falls to zero when the phase difference equals the value of n due to the symmetry of the ground 
state and the continuum state Yk ~ e tk r . For low harmonic orders, where the Coulomb effect is 
strong enough to modify the continuum state, is no longer a plane wave and the symmetry 
is broken. As a result, the phase different deviate from % and the nonzero ellipticity appears. 

To find in the frequency region where the Coulomb effect is unnegligible, we also compare 
our results with the analytical results in ref. [25] (also presented by the green curves in Fig. 
4(a)). For the harmonics in the range from the 21st to the 61st orders, the ellipticity of short 
quantum path obtained in our work is considerable and is in the same order of magnitude 
compared with that in (25). Note that in ref. (25), it is also in the same range from the 21st 
to the 61st harmonic order where their analytical result is lower than the numerical result. For 
harmonics higher than the 61st order, the ellipticity obtained in our work rapidly decreases to 
zero, and meanwhile the analytical results agrees well with the numerical results in ref. |25|. 
This comparison shows that, the Coulomb effect also contributes to the ellipticity for the low 
order harmonics, and neglecting the Coulomb effect would lead to the deviation of the result. 
For the long quantum path, our result is nearly negligible for all the harmonic orders compared 
with that in |25|. This is because the Coulomb effect is already very low for the long quantum 
path, where the continuum electron drifts far away from the core. 

To further investigate the role of Coulomb potential on the harmonic ellipticity, we also 
calculate the ellipticity with short-range potential Vj(r) = — \j\Jr 1 +aexp(— r 2 /s 2 ) [30], where 
s=1.6 and a=0. 102 to keep the ground-state energy unchanged. As shown in the inset of Fig. 4, 
the short-range potential approaches zero much closer to the core at about r=3 a.u. The obtained 
ellipticity with short-range potential is shown in Fig. 4(a) with dashed curves. Compared with 
the solid curves, the dashed curves are even sharper and drop to zero more quickly around the 
31th harmonic order. This is because the potential vanishes so rapidly that the trajectory of 
the associated continuum electron is already too far away for the short range of V s (r). In this 
relatively faraway area, the Coulomb effect is too weak to lead to the elliptically polarized HHG. 
However, for the HH of the lowest orders of short quantum path, the motions of the continuum 
electrons are confined in a very small range. Take the 21st harmonic order for example, the 
maximum distance of the electron trajectory from the core is only about 2.5 a.u. In this area, 
the difference of the Coulomb effects on these trajectories between the short- and long-range 
potentials is small. Therefore, the harmonic ellipticity of these orders is still high and become 
close to that obtained with soft-core potential. Similarly, as shown by the dashed curves in Fig. 
4(b), the phase differences also converge to the value of % much more rapidly due to the short 
range of the Coulomb potential. 



4. Conclusion 

In this paper, the role of the Coulomb potential on the generation of elliptically polarized HH 
from atoms driven by elliptically polarized lasers has been investigated analytically. It is found 
that the Coulomb effect makes a contribution to the harmonic ellipticity for low harmonic or- 
ders and short quantum path. We in detail analyze the influence of the Coulomb potential on 
the dynamics of the continuum state and find that the obtained ellipticity in our simulation 
originates from the break of symmetry of the continuum state. Moreover, the conclusions are 
further confirmed by comparing the results obtained with the soft-core potential and the short- 
range potential. We would like to further mention that this kind of Coulomb potential induced 
ellipticity is general for different species of core. Therefore, this mechanism can also, at least 
in part, result in the harmonic ellipticity in molecular HHG. In addition, with the analytical 
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Fig. 4. (a) The harmonic ellipticity and (b) the phase difference between the two orthogonal 
components as functions of the harmonic order. Inset: Comparison between the soft-core 
potential V(r) and the short-range potential V s (r). The analytical results in |25) are also 
presented in panel (a) with the green curves. 



model in SF-EVA, we have not only proposed the mechanism of the Coulomb effect on the har- 
monic ellipticity but also depicted the electronic dynamics of the Coulomb continuum states, 
which would be helpful for many other different topics in strong-field phenomena, such as finer 
molecular orbital tomography and ultrafast dynamic tracing. 
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